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Abstract. The transient dynamics of a quantum linear amplifier during the transition from damping to 
amplification regime is studied. The master equation for the quantized mode of the field is solved, and the 
solution is used to describe the statistics of the output field. The conditions under which a nonclassical 
input field may retain nonclassical features at the output of the amplifier are analyzed and compared to 
the results of earlier theories. As an application we give a dynamical description of the departure of the 
system from thermal equilibrium. 
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1 Introduction 

The master equation describing linear amplification or 
gain has been a prototype for discussing many questions 
in Quantum Optics. It was derived by the elimination of 
an unobserved environment using what has been termed a 
Born-Markov approximation . Thus its properties were 
mainly determined by physical considerations, but it ar- 
rived at a form later to be shown to be the consistent 
generator of dissipative time evolution in quantum the- 
ory; thus it is of the Lindblad form |2 . 

The master equation for linear amplification was ear- 
lier represented as the generic model for an optical ampli- 
fier or attenuator jSj. It also describes the onset of laser 
oscillations until the time when nonlinear saturation starts 
to affect the behavior. In the trapped ion context, when 
the ion trap potential is regarded in a harmonic approx- 
imation, the cooling by lasers may be considered as an 
attenuation described by the same equation j^- 

Its advantage is that it is exactly solvable, which allows 
us to follow the onset of gain or the damped approach to a 
steady state. The exact solution also allows us to evaluate 
the noise properties exactly and investigate the fading of 
nonclassical features of the initial state. 

In all applications so far, the amplifying and attenu- 
ating coefficients of the equation have been regarded as 
constants. This corresponds to the assumption that the 
population inversion is instantaneously reached and the 
evolution starts from an initial state experiencing no pre- 
vious evolution. Such theory, however, does not describe 
the transient dynamics of the linear amplifier, i.e. the dy- 
namics when the pumping field is switched on or off. In 
this paper we generalize the previous theory of linear am- 



plifiers to the case in which a smooth onset of amplification 
or attenuation takes place, and hence the amplifying and 
attenuation coefficients are time dependent. We solve the 
master equation with time dependent coefficients in terms 
of the characteristic function |5I6| and we use the solution 
to describe the transient dynamics of the linear quantum 
amplifier, one of the most widely used and common de- 
vices in quantum optics. 

We choose to consider a situation where the gain medium 
is switched smoothly from an attenuating regime to an 
amplifying one. This takes place within a time interval 
centered at some definite time, before which we have a 
damped situation. Thus the initial gain, which is normal- 
ized to unity of course, decreases first until the amplifi- 
cation coefficient changes sign and gain starts to grow. 
The model allows us to follow the solution through this 
point, and we can see how the time dependence affects 
the noise properties and the possibility to retain initially 
imposed nonclassical features of the system state. As we 
may expect, the situation is more complicated than the 
simple constant coefficient case. It is, however, possible to 
retain earlier results on amplifier added noise and quan- 
tum cloning limits by considering the appropriate limiting 
cases. 

Theoretical works on optical transients of physical phe- 
nomena which cause amplification of light have received 
in the past a huge deal of attention [7| . To the best of the 
authors' knowledge, however, this is the first analytic de- 
scription of the transient of phase-insensitive quantum lin- 
ear amplifiers. Therefore the results presented in this pa- 
per give a clear contribution to the fundamental research 
in the theory of lasers and optical amplifiers since every 
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linear amplifier or laser undergoes a transient behavior 
before stabilizing. 

Our results can be directly applied to describe tran- 
sients in technological applications based on linear ampli- 
fiers. Optical linear amplifiers are essential components of 
state-of-the-art optical networks. In order to attain best 
performances of the networks, however, it is crucial to an- 
alyze the behavior of linear amplifiers during power tran- 
sients causing fast switching on and off of the amplifiers 
|H] . Although the linear amplifiers currently used in optical 
networks do not need to operate at the quantum level, the 
recent development of quantum technologies such as quan- 
tum communication, quantum cryptography and quantum 
computation, is based on the implementation of networks 
containing nanodevices and nanocomponents operating at 
the quantum level. 

Recently a scheme for optimal cloning of coherent states 
with phase-insensitive linear amplifiers and beams split- 
ters has been proposed 9; . Recent advances in the field of 
nanoelectromechanical systems have paved the way to the 
realization of experiments close to achieve the quantum 
limited detection and amplification. In , position reso- 
lution very close to the quantum limit is obtained, demon- 
strating the near-ideal performance of a single-electron 
transistor as a linear amplifier. Another recent application 
of linear amplifiers consists in a method for reconstruct- 
ing a multimode entangled state Also quasiprobabil- 
ity functions have been shown to be measurable via direct 
photodetection of a linear amplified field jl2j. 

The paper is structured as follows. Section II sum- 
marizes the results of earlier investigations for easy com- 
parison with the present results. Section III presents the 
solution for time dependent coefficients and discusses its 
main properties. In Sec. IV we discuss the possibilities to 
retain nonclassical features in the output of the amplifier 
relating our results to earlier calculations. In Sec. V we 
discuss the emergence of thermal features in the solution. 
Finally Sec. VI concludes the discussion of the work. 



2 Review on phase insensitive narrow band 
linear amplifiers 

The simplest standard amplifier configuration consists of 
an assembly of N two-level atoms, of which N2 arc ex- 
cited and A'^i are unexcited, interacting with a single-mode 
quantum field. It is assumed that the field frequency is res- 
onant with the atomic frequency and that the population 
of the two-level atoms is partly inverted, i.e. N2 > Ni. 
This is the standard model of a laser; its linear operation 
regime describes an amplifier p. In the standard descrip- 
tion of linear amplifiers it is assumed that N2 and Ni are 
maintained approximately constant in time by some pump 
and loss mechanism. 

Starting from a microscopic description of the interac- 
tion between the two-level atoms and the quantum mode, 
it is possible to derive the following master equation for 
the field mode in the interaction picture 0] 

dp C 



A 



dt 



\a) ap — 2apa^ 



■ pa^a] 



aa) p 



2a) pa + paa^ 



(1) 



with a and annihilation and creation operator of the 
quantum harmonic oscillator and 



A 



C = 



7^ 



(2) 
(3) 



where g is the coupling strength of the interaction between 
the two-level atoms and the mode of the field, — Ni/^ 
(with « = 1, 2) is the pumping rate into the atomic level i, 
and 7 is a rate of the same order of the atomic linewidth. 

A relevant quantity in the dynamics is the linear gain 
(or damping) factor, describing the linear growth (loss) of 
energy in the mode. 



W = A-C. 



(4) 



When W > Q, the master equation describes a linear am- 
plifier, when W < Q \t describes a linear absorber. The 
constant A gives the noise provided by the spontaneous 
emission; this term is present even if the mode energy is 
initially zero. 

The time evolution of the amplitude of the field is de- 
scribed by the equation: 



{a),^t = {a{t)) = Gi/2e-*"»*(a), 



(5) 



where (a)i„ — {a{t — 0)), ojq is the frequency of the radi- 
ation field and the gain G is defined as 



(6) 



The gain is greater than 1 for linear amplifiers and smaller 
than 1 for linear attenuators. The solution of the Fokker- 
Plank equation for the Glauber-Sudarshan P representa- 
tion of the density matrix (P function) can be used to cal- 
culate the transformation of any incoming Pin (a) function 
by the amplifier: 



fout(a) = / daQP{a,t\ao)PiniaQ) 



where 

P(a,f|ao) 



TT'm{t) 



exp 



"0 



m(t) 



(7) 



(8) 



is the amplifier transfer function ,13] . The time dependent 
width is given by 



m{t) = A [G{t) - 1] /W. 



(9) 



This quantity represents the average photon number of 
the spontaneous emission field T. Note that, for a linear 
amplifier, the gain grows asymptotically to infinity for t 
00, and so does the width [3]. For an absorber, on the other 
hand, the asymptotic value of the width m{t) is finite. 
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By using Eqs. Q and one can calculate the noise 
of the output field, defined as the symmetrically ordered 
fluctuations of the field mode, 14 : 

l^alLt = ^{a^a + aa^)oiit - {a)out{a^)out 



= G\Aa\l+m{t)--{G-l) 
^G{\Aa\l+A), 



where 



•^"2 [A-C 



1 



(10) 



(11) 



is the equivalent noise factor, or amplifier added noise, in- 
troduced by Caves ^Ij. This quantity describes the fluc- 
tuations of the internal modes of the amplifying medium. 
Since the input field and the internal modes of the am- 
plifying medium are uncorrelated, their fluctuations add 
in quadrature and they are both amplified. The minimum 
value of the added noise for infinite gain, e.g. for t oo, 
is given by the Caves limit: 



1 fA + C 



2\A~C 



1 

^2' 



(12) 



where the excess noise factor 9 gives the initial mean num- 
ber of excitations of the internal modes of the medium, 
and therefore approaches zero when the initial tempera- 
ture of the amplifying medium vanishes, T — > 0. 

It has been demonstrated that the output field of a 
phase insensitive narrow band linear amplifier may possess 
nonclassical features only if the input field is nonclassical. 
However, as the amplifier gain increases, any nonclassical 
feature of the light, which was present in the input field, 
tends to be lost. In particular, subPoissonian statistics 
and squeezing are lost when the gain G exceeds the value 

2 [mg. 

In the next section we present a theory describing 
the transient regime of the amplification process. In other 
words, we will drop the assumption that N2 and A^i are 
constant and we will describe the onset of the amplifica- 
tion process from an initial damping regime. Our aim is 
to study the transient dynamics and to investigate how 
the results for the standard amplifier, described in this 
section, are modified. 



3 Transient regime of linear amplification 



case, the time evolution of the field mode is described by 
a master equation of the same form of Eq. (Q, but with 
time dependent coefficients A{t) and C{t). Similarly, if one 
switches off the external pump, the ratio of the atomic 
populations N2/N1 will tend to the Boltzmann factor and 
the amplification process will eventually stop, the system 
approaching its thermal equilibrium. In the following we 
consider the first of these two situations, i.e. the onset of 
amplification due to the creation of population inversion 
in an initially damping medium. In more detail, we con- 
sider the case in which 

2o2 e{t-to) 

A(i) = -Vr2(<) = A^— rm+B, (13) 

where e is the rate of change of the pumping coefficients, 
that is the amplification onset rate. We assume that at 
t — —00 the state of the ensemble of two- level atoms con- 
stituting the amplifying medium is thermal, that is 



N2{t^-oo) A{t^ -00) 



B 



Ni{t 



-00) C{t 



A + B 



(15) 



-huiQ/ksT 



which implies B/A — (e 



_ I huJa/kBT 



— 1) = um, with um 



mean number of excitations of the medium. We assume 
that the state of the amplifying medium practically does 
not change in the time interval —00 < t < 0. 

Under these conditions the asymptotic gain factor W{t) 
takes the form 

W{t) = A(t) - C{t) = Atanh[e(i - to)]- (16) 

Note that, for t < t^, W{t) < Q and the system behaves 
as an absorber, while for t > to, W{t) > and the system 
behaves as an amplifier. Therefore Iq indicates the time at 
which the amplification process begins. From Eq. H16|l we 
infer that the constant A is the asymptotic gain factor. 
Finally we stress that, for t —t 00, 

N2 (t oo)/Ni{t ^ 00) = A{t 00) /C{t 00) 
= {A + B)/B = Ni{t^ -00)/ N2{t^ -00), (17) 

i.e., the population of the exited (ground) state tends 
asymptotically to the initial population of the ground (ex- 
cited) state. 



Previous work on linear amplifiers deals with a situa- 
tion in which the amplifying medium, e.g. an assembly 
of two-level atoms, satisfies the population inversion con- 
dition required to amplify an input field. In order to reach 
the inverted population condition it is necessary to pump 
the atoms from their initial thermal condition till the 
point in which N2 > Ni. During this transient regime 
the pumping rates to levels 2 and 1 (r2(t) and ri(t)) 
change with time till they reach a stationary value for 
which N2/N1 (X /ri > 1 (amplification regime). In this 



3.1 The master equation and its solution 

The master equation describing the transient behavior of 
the amplification process, in the interaction picture, is the 
following 



dp 
d7 



2 

C'ir) 



aa) p - 



2a) 



pa 4- paa^ 



\a) ap — 2apa^ 



(18) 
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with 



A'(t) = A' — 

^ ' p(t 



e 



e(T-ro) ^ g-(T-ro) 

g-(T-To) 



(r-To) ^ e-(T-To) 



(19) 
(20) 



In the previous equations we have introduced the relevant 
physical dimensionless parameters t = et, tq = eto, A' = 
A/e and B' = B/e. As we will see in the following, the 
parameter A' , which is the ratio between the asymptotic 
gain factor and the rate of onset of the amplification, plays 
a central role in the system dynamics. Indeed both the 
gain G(t) and the added noise depend crucially on this 
parameter. 

Following the method developed in [S] we solve the 
master equation given by Eq. (|18|l in terms of the quantum 
characteristic function (QCF) jl^, defined through the 
equation 



PS 



The solution reads as follows 



XriO 



Xo 



-i{ujo/e)T 



(21) 



(22) 



where xo is the QCF of the initial state of the field, luq is 
the field frequency and the G(t) is the gain, given by 



G(r) 



(23) 



The quantity /^(t). 
follows 



appearing in Eq. H22I) is defined as 



^{r) = \[G{t)] IjGir')]-' [G'(r') + A'(r')] dr' . 



(24) 




Fig. 1. Time evolution of G(r) for A' = 2, 3, 4, 5.5 (increasing 
values of A! correspond to decreasing thickness of the line) and 

To = 5. 

calculations, it turns out that the P function has the same 
form of Eq. but with G(t) given by Eq. (|2Sl and 
m(r) = [G(r)-l]/2 + Z\(r). 

Inserting Eq. (|16|l into Eq. H23|l and carrying out the 
integration yields 



G(r) 



cosh(T — To) 
cosh(ro) 



A' 



(26) 



It is worth underlining that the solution given by Eq. (|22() , 
with the help of Eqs. f^ - f^ . holds whatever the ex- 
plicit time dependence of the coefhcients A(t) and G(r), 
appearing in Eq. (|18|l . is. The case considered in the pa- 
per [Eqs. (|13|l - H14(l ] has been chosen to illustrate the tran- 
sient dynamics in a physically reasonable and well justi- 
fied model. Indeed Eqs. (|13|I - H14() describe a situation in 
which from an initial condition in which A^i > A^2, the 
populations of the two-level systems constituting the am- 
plifying medium pass smoothly to the inversion condition 
.^2 > -^1 necessary for amplification. In passing, we note 
that the hyperbolic tangent time dependence is one of the 
most commonly adopted phenomenological models in the 
description of transients of physical systems. 

Starting from Eq. (|22|l . one can calculate the Wigner 
function, the Glauber- Sudarshan P function, and the Husimi 
Q function by means of the relation jl5| 



It is not difficult to prove that, for tq — 0, and in the limit 
of infinitely fast onset of the amplification (e cxd), the 
gain function tends to G(t) = e'^* = e^* [see Eq. Q]. 

In Fig. n we plot the gain G(t) for four increasing val- 
ues of A! . As clearly shown in the figure for increasing 
values of A! the values of the gain in proximity of the 
amplification time tq become smaller and smaller. This is 
because small values of A! correspond to small values of 
the asymptotic gain factor or, equivalently, to a very slow 
amplification onset rate. In general, the gain decreases for 
times smaller than tq and, as expected, starts to increase 
after the amplification sets in, even if G(t) > 1 only for 
times r > 2ro. Note that the standard theory of linear am- 
plification predicts that for a linear amplifier it is always 
G(t) > 1 (see Sec.|21l. However, if one takes into account 
the transient regime characterizing the initial dynamics 
of every linear amplifier it turns out that there exist an 
interval of time at which, although > 0, there is still 
no gain. The reason why the gain becomes greater than 1 
only after the time 2to is that for < r < tq the system 
is in a damping regime, and hence the gain decreases. It 
takes exactly another interval of time tq after the onset 
of the amplification process to undo the initial decrease 
in the gain. At r = 2to we have G(2to) = G(0) = 1 af- 
ter which the gain increases monotonically. This is a new 
feature brought to light by our theory. 

Let us focus on the quantity A{t). Inserting Eq. (|26|l 
into Eq. if^ we get 



d2ex.(0exp(ar-a*0e(fl«l'/2), 



(25) 

where -p = —1,0,1 corresponds to the Q, Wigner, and 
P functions, respectively. In particular, carrying out the 



A{t) = G(t) 



A' -t- 2B' 



C0Sh(T' — Tq) 



cosh(To) 



-A' 



dr' 



G{t) Ia'{t~tq). 



(27) 
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For tq — and for each A' real, with A' ^ 1, we have 
sinh r 1 



lA'iT) = 



(28) 



(A' - 1) (coshr)^'"! 

1 - A'/2,3/2 - A'/2- (coshr)^], 

with l-A'/2, 3/2-A72; (coshr)^] being the hyperge- 
ometric function of the variable x — (coshr)^. For A' — 1 
the integral appearing in Eq. (|27|l is simply equal to 



/i(t) = 2arctan(e^). 



(29) 



The mathematical expression of the added noise in the 
special case of integer values of A' is discussed in Appendix 
A. For To 0, one gets 

Z\(r) = G(r) ^' + ^^' cosh(To)^' 

[/A'(r-ro)+/A'(To)], (30) 

where Ia'{t — tq) and Ia'{to) are obtained from Eq. H28|l 
by substituting for the variable r the expressions r — tq 
and To, respectively. 



3.2 Noise of the output field 

From the QCF solution given by Eq. H22|l we can easily 
calculate the mean values of observables of interest, e.g. 
those characterizing the output field statistics, by means 
of the relation ^HJ 



(aT™a") 



d 



x(0 



(31) 



We look first of all at the symmetrically ordered fluctua- 
tion, as defined by Caves 14 



\Aa\l^,^G{T)\Aa\l + 
= G(r) [\Aa\l H 

The added noise is given by 

A = A{t)/G{t), 



Air) 



(32) 



(33) 



with A{t) given by Eq. (jSOJ. This quantity is clearly dif- 
ferent from the one given by Eq. for the standard 
linear amplifier case. It is possible to show that for infi- 
nite gain, e.g. for t oo, this quantity is always greater 
or equal to Ac^ Ac being the Caves limit for an infinitely 
fast onset of the amplification process. 

In order to derive the Caves limit from Eq. 133|l . we 
note that for r — > oo, Ia', as given by Eq. (|28|l . tends to 

CHI 



TOO 

^A' — 



R r{A'/2) 



2 r[(A' + l)/2]' 
Therefore the asymptotic value of the added noise is 



(34) 



A = 



A' + 2B' 



Ia' 



A' + 2B' r{A'/2) 



2 r[{A' + l)/2] 



r{A' /2 + 1) 

rp' + l)/2]- 



(35) 
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-0.4 
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Fig. 2. Damping and amplification of an initial input field 
having (a(t = 0))=10, for tq = 4, A' = 0.5 and B' = 0.5 ■ 10"^. 
The dashed lines indicate signal width Z\(r)^'''^. The insert 
shows the gain factor VK(r) = Air') — C{t) in the same time 
interval and for the same values of the parameters. 



The Caves limit is obtained for an infinitely fast onset of 
the amplification process, that is for e —^ oo, i.e. A' — > 0. 
Substituting A' = into Eq. H35() . and remembering that 
r(l) 1 and r(l/2) = Vtt, one gets 



1 B 1 



(36) 



that is the Caves limit [see Eq. (|12|l ]. From Eq. H35I) one 
can see that, for fixed values of the initial mean number 
of excitations of the medium nm — B/A, Ac is actually 
the smallest asymptotic value of the added noise. 

A careful analysis of the noise at the output field, as 
given by Eq. 1)32(1 . shows that this quantity, as one would 
expect, increases monotonically with time whatever the 
initial state is. In Fig. [5] we show the dynamics of an input 
field which is damped for r < tq = 4 and then amplified 
for r > To = 4. The figure shows that the width of the sig- 
nal always increases. In the following section we will study 
in more detail the transient dynamics for different types 
of input fields. Since we are dealing with phase insensitive 
amplifiers/ absorbers, it turns out that it is not possible 
to generate nonclassical states from classical input fields. 
However, it is possible to analyze how the conditions to 
retain initial nonclassical features are modified due to the 
transient dynamics. In addition we will look at the field 
statistics by explicitly calculating the time evolution of 
the Wigner function. 



4 Nonclassical properties of the output field 

4.1 Squeezing and subPoissonian statistics 

Let us begin studying the conditions for which the output 
field can retain squeezing properties when the input state 
is a squeezed state of the electromagnetic field. We define 
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the dimensionless quadratures of the field as foUows 

1 



V2 



(a - at) . 



(37) 
(38) 



The squeezed states satisfy the minimum uncertainty re- 
lation AuAv = 1/2, but are characterized by an unequal 
distribution of the quantum fluctuations 



Au = 



Av 



1 



(39) 



with s 7^ 1. Introducing the rotating coordinates 



u 

V 



u cos{uj^t) — V sin(cL;o^), 
V cos(aJo^) + sm{LOf)t), 



and using Eq. H31(l we get 



(40) 
(41) 

(42) 
(43) 

with A given by Eqs. (|33|l and (|27|1 . For an input squeezed 
state having s < 1, the output can remain squeezed if and 
only if 

G{t)[{Au)1,^A]<\. (44) 



{Au)1^,=G{t) [{Au)l,+A], 
{Ai)l^,^G{T) [{Ai)l, + A], 



It is possible to show that the maximum allowed value of 
the gain G(t), in order to retain squeezing at the output, 
decreases with A' and B' . In other words, for increasing 
values of A' and J5', one can retain squeezing only for 
smaller and smaller values of the gain (less efficient am- 
plification). Having in mind Eq. (|30|l . one finds that, for 
To = 0, the output field is still squeezed if the gain satisfies 
the following inequality 



G{r) < — 



{A' + 2B')Ia'{t) 



(45) 



It is worth recalling the standard result for phase insen- 
sitive linear amplifiers, which states that the upper limit 
for the gain compatible with squeezing at of the output 
field is G = 2, the magic number for photon cloning TS . 
The analysis of the behavior of the gain in our case is 
more complicated, since the r.h.s. of the inequality H45|l 
depends on time. A numerical study shows that, although 
for certain time intervals, the r.h.s of the inequality may 
be greater than 2, in these time intervals G(r) is always 
smaller than 2. Hence, also in the case studied in this 
paper G = 2 constitutes an upper limit for retaining non- 
classical features in the output field. 

Let us now look at the dynamics of an input Fock 
state. We recall that one of the nonclassical features of 
such states is that their statistics is subPoissonian. Simi- 
larly to what we have done for the squeezed states we an- 
alyze the requirements to retain subPoissonian statistics 
at the output field. To this aim we introduce the Mandel 
parameter Q P 



Q 



in) 



1. 



(46) 
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6 TQ 



(b) 
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G(tq)=1.32 
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0.6 IQ 0.8 



1 



Fig. 3. Mandel parameters of the output fields for an initial 
Fock state |no = 5) in the cases A' — 0.05 (a), and A' = 1 (b). 
We have set ro = and B/A= 10"^ in both (a) and (b). We 
indicate with tq the instant of time at which Q = 0. In the 
upper-left corner we indicate the corresponding value of the 
gain at r = tq: G{tq). 



This quantity gives an indication of the statistics of a 
quantized field. For a Fock state Q takes its lowest value 
Q = — 1 while for a coherent state Q is equal to . There- 
fore, values of Q < indicate subPoissonian statistics, 
while Q = indicates Poissonian statistics and Q > 
superPoissonian statistics. Using Eq. lp?T|l we derive the 
time evolution of the Mandel parameter as follows 



Qo 



[G(r)]»„JQ,„-(n) 
{n)oiit 



(47) 



where {n)out and are the mean number of photons 

of the output and input fields, respectively, and 



{n)out = G(r)(n)„, + i [G(r) - 1] + A{t) 
= G{T){n),r,+m{T). 



(48) 



For a Fock input state Qm = — 1 and {n)in — hq, hence 
the condition for having subPoissonian statistics at the 
output is 



(n) out {{n) out + 1) < G^{T)nQ{nQ + 1). 



(49) 
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A numerical analysis shows that, as for the squeezing, also 
for the subPoissonian statistics, we obtain the limit 2 for 
the gain typical of the standard theory of linear amplifiers. 

As an example, in Fig. El we compare the Mandel pa- 
rameters of the output fields for an initial Fock state |no = 
5) in the cases A' = 0.05 (fast onset of the amplification 
and/or small value of the asymptotic gain) and A' — 1. 
From the figure one sees that, for A' = 0.05, one can re- 
tain subPoissonian statistics of the output field for higher 
values of the gain compared to the A' = 1 case. 



4.2 Wigner function 

Let us now have a look at the complete statistic of the 
output field, by means of the Wigner function. Inserting 
Eq. (|^ into Eq. , and putting p — we get 



(50) 



Inserting the inverse Fourier transform of Eq. 125(1 into 
Eq.ljHOI) gives 



Wr(a) 



rf^aoW^o(ao)- 



exp 



"o)4'-b'(T,Q,ao)C 
&(r,a,ao) |' 



A(r) 



A{t) 



d?aoWr{a\aoyVo{ao), 



with 



6(t, a, ao) = a - aoG^/'^{T)e''^°^'' . 



(51) 



(52) 



In the derivation of Eq. H51|) we have used the property 
that the Fourier transform of a Gaussian is a Gaussian. 
The quantity VFT-(a|Q;o) is the propagator which, for r — > 
0, tends to the delta function 5(a — ao). 

If the state of the input field is a coherent state |ao), 
than the Wigner function of the output state reads as 
follows 



exp 



A{r) + l/2 



A{t) + 1/2 



(53) 



The Wigner function of the output state is therefore a 
Gaussian. Having in mind the time evolution of G(r) [see 
Ea. H26l) and Fig^, one realizes that, in a frame rotating 
with the frequency loq, the Wigner function of an input 
coherent state |ao)j with ag ^ 0, moves towards the center 
of the phase space for r < tq and then moves away for 
r > To, while its width continuously increases. 



(a) 



T=0 



x=6 
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2.5 



-2.5 
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1=12 



x=18 




-7.5-5 -2.5 2.5 5 7.5 



-7.5-5 -2.5 2.5 5 7.5 
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Fig. 4. Contour plots of the Wigner function at different times 
r for an input squeezed state with r = 1. We have put tq = 4, 
A' — 0.1, ub = B/A = 0.1 (a). Variance of the quadrature u 
as a function of time, for the same values of the parameters 
(b). 



We now consider the case of an initially squeezed state. 
The initial QCF for squeezed coherent state is 



XoiO = exp 



-iea-re-^^5,|2 + i(ras + eao) 



Here Cr — cosh(r) and Sr — sinh(r) , ao is the displace- 
ment of the input field and z = re~"^ is the squeezing 
argument. 

For an input squeezed vacuum state (ap = 0), with 
squeezing angle = 0, the Wigner function at time r 
takes the form 



1 r°° 

7^ J~oo 



exp 



(55) 
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This Fourier transformation can be calculated with the 
method used in ^17., the result being 



Wria) = Mexp 



~2ai 



exp 



2Z\(t)G(t)(C2,. -^2r)-l 



(56) 



Here, ax and ay are the real and imaginary parts of a, and 
M is a time dependent normalization constant. We note 
that this result is consistent with the Eqs. (|42(I - H43() used in 
Sec. l4.1l to study the time evolution of the quadratures of 
the field for an initial input squeezed state. Contour plots 
showing the time evolution of the Wigner function are 
shown in Fig.^(a), while in Fig. 0](b) we show the time 
evolution of the squeezing of the quadrature amplitude u. 



5 Departure from thermal equilibrium 

The analytic approach we have described in the previous 
section to analyze the onset of the amplification process, 
can be used to study how a system departs from an initial 
thermal equilibrium situation. In more detail, we consider 
the case in which the medium, modelled as an ensem- 
ble of two-level atoms, is initially in thermal equilibrium 
with the mode of the quantized field. The ratio N2/N1 
between the populations of the excited and ground states, 
respectively, is therefore given by the Boltzmann factor 
N2/N1 = e-^'^o/'^B^. The state of the field is a thermal 
state at T temperature. 

At T = one switches on pumping lasers which change 
the population of the two-level atoms until the condition 
of population inversion, necessary for the onset of the am- 
plification process, is reached. The pumping lasers alter 
the initial condition of equilibrium between the medium 
and the system (the field mode). In order to study how the 
system departs from the condition of thermal equilibrium 
with the two-level atoms medium, we use the solution of 
the Master Equation (|18|l to calculate the time evolution 
of the Wigner function of the field. For an initial thermal 
state, the QCF at time r has the form 



e-^MI«l exp 



[n)out + l]G{T)\e 



(57) 



Inserting this equation into Eq. H25|l one gets the following 
expression for the Wigner function at time t 



1 



Wr{a) = - 



1 



TT {n)out + 1/2 



cxp 



{n)out + l/'2 



(58) 



In the last two equations, {n)out is the number of pho- 
tons of the output field, as given by Eq. The Wigner 
function of Eq. (|58l) is the Wigner function of a thermal 
state at a temperature T{t) which varies with time. The 
medium is not in thermal equilibrium anymore, since the 
pumping lasers change the two-level atoms population un- 
til the population inversion condition is reached. However, 
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Fig. 5. Time evolution of the temperature of the system for 
A' = 1 (dotted line) and A' = 0.05 (solid line). For both graph- 
ics we have set ojq = W^'^Hz and ub = B/A — 10^. The box 
in the top left corner is the ratio between the populations of 
the two-level atoms. This quantity does not depend on A' and 
B' separately, but on the ratio ns = B/A only [see Eqs. I|13|l 
and llT^ I. 



Ea. (|58(l shows that the medium plus the pumping lasers 
behave, as far as the system (mode field) is concerned, as a 
thermal reservoir at varying temperature T'(t), as known 
from the theory of laser cooling 20 . Having in mind that 
2(n)oMt + l = coth [fiuJo/kBT{T)] [Ij, and using the relation 
arcoth(x) — [lii{x + 1) — ln(x — l)]/2 (for x"^ > 1), we can 
express the time evolution of the temperature as follows: 



T(r) 



hujQ 



{ln[(n)o„t + 1] -ln[(n)o„t]} V 



(59) 



Figure El compares the behavior of T(r) for the two cases 
A' = 1 and A' = 0.05. The figure shows clearly that in 
both cases the temperature of the system is constant dur- 
ing the damping regime and it starts to increase when 
approaching the population inversion at r = tq, i.e. when 
the change in the population ratio N2/N1 becomes con- 
siderable. The increase in the temperature is much higher 
for higher values of A' , since in this case the asymptotic 
gain is also higher. 

In order to characterize further the departure from 
the initial condition of the system we look at the von- 
Neumann entropy of the field. We use the result obtained 
by Agarwal to calculate the dynamical entropy for a 
state of the form given by Eq. lf55|) 



Sir) = ks {[{n)out + l]ln[(n)o„f - 
- (n)o„t In [(n)o„t]} . 



(60) 



From direct inspection in the previous equation one sees 
that, similarly to the dynamics of the temperature, the 
entropy remains approximately constant for r < tq and 
begins to increase when r ~ tq. In Fig. |^ we show how the 
entropy increase rate (in units of fc^) changes with A'; in 
the box in the top left corner the dynamics of the entropy 
for three example values of A' is shown. For increasing 
values of A' , the linear increase in the entropy due to the 
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laser pumps up the atoms of the medium till the condi- 
tion of population inversion is reached. This is an example 
of dynamic departure from a thermal equilibrium condi- 
tion that can be studied analytically. We analyze the time 
evolution of the temperature and of the von-Neumann en- 
tropy on the characteristic parameters of the linear am- 
plifier. We find that, as known from the theory of laser 
cooling, the medium plus the pumping lasers behave, as 
far as the system is concerned, as a thermal reservoir at 
varying temperature. We find that the entropy increase 
rate depends crucially on the asymptotic gain. 



A' 7 Acknowledgements 



Fig. 6. Dependence of the entropy increase rate on 
A'. The entropy increase rate is defined as AS/Arks = 
[S(r = 14) - S{t = 10)] /4fc_g. In the box in the top left corner 
the time evolution of the entropy for the three exemplary val- 
ues A' — 1, A' — 0.5, and A' — 0.05 is shown. In all the plots 
we have set tq — 8, and ub ~ B/A = 10. 



amplification process is faster and faster. This result is in 
accordance with the behavior of the temperature of the 
system. In fact, smaller values of A' correspond to smaller 
asymptotic gain and therefore less efficient amplification 
processes. 



6 Conclusions 



In this paper we have discussed the dynamics of a quan- 
tum linear amplifier during the onset of the amplification 
process. For an amplifying medium consisting of an assem- 
bly of two-level atoms, our theory describes the dynamics 
of the output field when the medium passes from a con- 
dition in which the population of the atoms is thermal, 
to a condition of population inversion characterizing the 
amplifying regime. 

We have solved exactly the master equation describing 
the transient dynamics of the linear amplifier in terms of 
the quantum characteristic function. The solution is used 
to investigate conditions under which an input nonclassi- 
cal field may retain nonclassical features at the output of 
the linear amplifier. We derive the analytic expressions for 
the output noise, as well as for the squeezing, the Mandel 
parameter and the Wigner function of the output field, 
and we use them to characterize completely the transient 
dynamics of the output field. 

Our results are compared with earlier theories of phase 
insensitive linear amplifiers which rely on the assumption 
that the population inversion is instantaneously reached, 
i.e. neglecting the transient regime. We show that also 
for a slow onset of amplification, the gain G{t) has to be 
smaller than 2 (the cloning magic number) in order for 
the output field to retain initial nonclassical properties. 

We conclude the paper analyzing the situation in which where 
the initial mode of the field and the two- level atoms medium 
are in thermal equilibrium at T temperature. An external 
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Appendix A 

For integer values of A' the integral Ia' defined in Eq. H27|l 
is given by [T^ 



sinh T 



1 



2771 — 1 (coshr)^ 



fe=i 



I'. 



2m+l 



r{m - k)r{m - 1/2) 
sinh r 1 



(61) 



2m (coshr) 



2m 



1 I V ^(-"-fe)^('" + ^/^) (co.hr)^^ 
+ r{m)r{m - k + 1/2) ' 

(2m- 1)" 



(2m)! 



■ arctan(sinh t) 



(62) 



In this appendix we show that the two equations written 
above are special cases of Eq. (|^ . 

For A' = 2m > 1 an even integer, the hypergeometric 
function reduces to a polynomial of order 1 — A' /2 



F[-m,b,c;z] = ^ 

k=0 



{-m)k{h)k 



(c)a 



fc!' 



r{z + fc) 
r{z) ' 



(2)0 = 1. 



(63) 



(64) 
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Using Eq. H63|l and the following properties 



14, 
15, 



CM. Caves, Phys. Rev. D 23, 1817 (1982) 

S.M. Barnett and P.M. Radmore, Methods in Theoretical 

Quantum Optics (Clarendon Press, Oxford 1997) 

M. Abramowitz and I. Stegun, Handbook of Mathematical 

Functions (Dover, New York 1965) 

K. Matsuo, Phys. Rev. A 47, 3337 (1993) 

I.S. Gradshtein and I.M. Ryzhik, Tables of Integrals, Series 

and Products (Academic Press Inc., San Diego 1994) 

G.S. Agarwal, Phys. Rev A 3, 828 (1971) 

S. Stenholm, J. Opt. Soc. Am. B 2, 1743 (1985) 



r{z + i) = zr[z) 




16, 



equation (|28|l reduces to Ea. ljtilll . 



17. 
18. 



For A' = (2m + 1), we obtain Eq. from Eq. ||2Hl 
by using the properties 



19, 
20, 



- i sinh(T)F[l, 1 - A72, 3/2 - A! jl; (coshT)^] 
= i^[l/2 - A! 12, l/2,i/2 - A! 12; (coshr)^] 
= i^[-m, 1/2, 1 - to; (coshr)^], 



and 



F[-TO, 1/2, 1 - to; (coshr)^] 

= ^r{\ - TO)(-coshT)™P;^'(zsinhT), 



(65) 



where 




(66) 



with Pm(z) Legendre polinomials. 
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